Generalized Fluctuation Theory based on 
the reparametrization invariance of the microcanonical ensemble 
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Abstract 

The main interest of the present work is the generalization of the Boltzmann-Gibbs distributions and the fluctuation theory 
based on the consideration of the reparametrization invariance of the microcanonical ensemble. This approach allows a novel 
interpretation of some anomalous phenomena observed in the non extensive systems like the existence of the negative specific 
heats as well as possibilities the enhancing of some Monte Carlo methods based on the Statistical Mechanics. 



I. INTRODUCTION 



We have been witness in the last years of an important 
reconsideration of the Foundations of the Thermodynam- 
ics and the Statistical Mechanics originated from the in- 
terest to use this kind of description in systems which are 
outside the context of their traditional applications. An 
important area of the nowadays developments is found 
in the Thermo-statistical description of the called non 
extensive systems. 

The non extensive systems are those non necessarily 
composed by a huge number of constituents, they could 
be mesoscopic or even small systems, where the char- 
acteristic radio of some underlaying interaction is com- 
parable or larger than the characteristic linear dimen- 
sion of the system, particularity leading to non existence 
of statistical independence due to the presence of long- 
range correlations. In the last years a significant volume 
of evidences of thermodynamic anomalous behaviors are 
been found in the study of the non screened plasmas and 
the turbulent diffusion, astrophysical systems, nuclear, 
molecular and the atomic clusters, granular matter and 
complex systems 0, i g i, i, 1, 1, 1, E3, E3 • 

It is usual to find in this context equilibrium thermody- 
namic states characterized by exhibiting a negative heat 
capacity, that is, thermodynamic states where an incre- 
ment of the total energy leads to a decreasing of the tem- 
perature. This behavior is associated to the convexity 
of the entropy and the corresponding ensemble inequiv- 
alence and can appear as a consequence of the small 
character of the systems (during the first-order phase 
transitions in molecular and atomic or nuclear clusters 
[13, EE EH) or because of the long-range character of the 
interactions (see in the astrophysical systems [HI). 

The existence of such anomalous thermodynamical be- 
haviors demands new developments within the Equilib- 
rium Statistical Mechanics and the Thermodynamics. 
An important contribution to the understanding of such 
phenomena is found in the formulation of the Micro- 
canonical Thermo- statistics by D. FT. E. Gross [7|, a the- 
ory which returns to the original basis of the Statisti- 
cal Mechanics with the reconsideration of the well-known 



Boltzmann epitaph: 



S = logVF, 



(1) 



and allows the study of phase transitions in systems out- 
side the thermodynamic limit. Since W is the number of 
microscopic states compatible with a given macroscopic 
state, the Boltzmann entropy |1| is a measure of the size 
of the microcanonical ensemble. Within Classical Statis- 
tical Mechanics W is just the microcanonical accessible 
phase space volume, that is, a geometric quantity. This 
explains why the entropy (Q]) does not satisfy the con- 
cavity and the extensivity properties, neither demands 
the imposition of the thermodynamic limit or a proba- 
bilistic interpretation like the Shannon-Boltzmann Gibbs 
extensive entropy: 



(2) 



Such geometric character of the Boltzmann entropy leads 
naturally to the development of geometric formulations of 
any statistical formulation derived from microcanonical 
basis. 

Motived by the successes of the Gross' formulation 0], 
we have performed in the refs.[l7iEU an analysis about 
the geometric features of the microcanonical ensemble. 
Particularly, we have shown in the previous paper [ItJ 
that the microcanonical description is characterized by 
the presence of an internal symmetry whose existence is 
related to the dynamical origin of this ensemble, which 
has been coined by us as the reparametrization invari- 
ance. Such symmetry leads naturally to the develop- 
ment of a non Ricmannian geometric formulation within 
the microcanonical description, which leads to an unex- 
pected generalization of the Gibbs canonical ensemble 
and the classical fluctuation theory for the open systems, 
the improvement of Monte Carlo simulations based on 
the canonical ensemble, as well as a reconsideration of 
any classification scheme of the phase transitions based 
on the concavity of the microcanonical entropy. 

I shall continuous in the present paper the study of 
the geometric aspects related to the existence of the 
reparametrization invariance. The present work ad- 
dresses the study of those equilibrium situations char- 
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acterized by the consideration of several control param- 
eters. The main interest will be focussed in the general- 
ization of the well-known Boltzmann-Gibbs distributions 
and the fluctuation theory related with the consideration 
of the reparametrization invariance ideas. I realize while 
writing this work that some of the present ideas have been 
also rediscovered by Toral in ref. [19j | by starting from the 
Information Theory. 



II. ISOLATED SYSTEM 

Let us consider an isolated dynamical system whose 
macroscopic description could be carried out by starting 

from microcanonical basis. Let I = jii,/^--- In\ be 
the set of all those fundamental physical quantities de- 
termining the microcanonical distribution function ui m : 



UJ, 



(I) 



1 



n(i) 



i - 1 



(3) 



which are simultaneously mensurable, that is p, I 

0, being fl (I) = Sp |(5 I — I j the microcanonical par- 
tition function. Let us suppose that the eigenvalues of 
the set / belong to certain subset IZi of the Euclidean 
n-dimensional space lZ n . 

Let us consider another subset TZq C lZ n which is dif- 
feomorphic equivalent to the subset IZi, that is, there 
exists a diffcomorphic bijective map : IZi — > TZq with 
a non vanishing Jacobian: 



r<90 fe (I) 



det 



dP 



(4) 



The above identity ((6|) expresses an internal symme- 
try which will be referred as a reparametrization invari- 
ance of the microcanonical ensemble. The reparametriza- 
tion invariance of the microcanonical distribution func- 
tion ui m leads to the reparametrization invariance of 
the expectation value of any microscopic quantity A, 

A = Sp^Alu„^ =>■ A{I) = A(Q), that is, the expec- 
tation values of the microscopic observables behaves as 
scalar functions under any reparametrization change 0. 

Mathematically speaking, it is said that the Euclidean 
subsets IZi and TZq are identical representations of cer- 
tain abstract space 3. The totality of the reparametriza- 
tion changes among the admissible representations of the 
space 5 constitutes a group, the Group of Diffeomor- 
phisms of the space 3. The reparametrization invariance 
of the microcanonical ensemble means that such descrip- 
tion provides the same macroscopic Physics for all ad- 
missible representations of the abstract space 3. 

The microcanonical partition function f2 defines an in- 
variant measure d[i = Q (0) dO = (I) dl which pro- 
vides the number of microstates W a = J dfj, belonging 
to a small subset 7r CT of certain coarsed grained partition 



V = 



IK = 9 



(7) 



of the abstract space 9, and allows us to introduce the 
Boltzmann entropy S*b [tt ct ] = lnlTer. Hereafter, I will 
consider that the number of macrostates W a grows ex- 
ponentially with the increasing of the number n of the 
degrees of freedom of the interest system, so that the 
thermodynamic limit: 



lim 



(8) 



for every point I G TZi- It is said that the map tp de- 
fines a reparametrization change of the Euclidean subset 
IZi. The set of functionals derived from the above map 



e, © = I© 1 ,© 2 ...©"} 







, determines the same 



microscopic Physics of the set /: every eigenstate 



of J, I\* a ) = I a 



©, ei* Q ) 







^q), is also an eigenstate of the set 

|*a> = e[J a ]|* a ) = 9 a |* a ), and 



the physical quantities fc are also simultaneously men- 



surable, 



Qk Qm 



0. There is nothing strange that the 



set determines also the same macroscopic Physics of 
the set /. Taking into account the well-known property 
of the Dirac delta function: 







<90 


-l 




0-0 




5 


i- i 






dl 







fi(0) = 



90 



81 



and therefore: 



fi(0) 



0-0 



Q(I) 



I - I 



n(j), 

(5) 
(6) 



exists and it is independent of the nature of the coarsed 
grained partition ir a of the space 3 whenever every sub- 
set TT a be small. This hypothesis allows us to take the 
reduced entropy s as a continuous scalar function de- 
fined on the space 3. When the interest system is large 
enough, the number of microstates W can be estimated 
in the practice by W\ ~ f2 (/) 6Iq, where 5Iq is a small 
constant volume. The estimate W\ differs from the one 
obtained in other representation Wa — ^ (0) S@q. How- 
ever these estimations should lead to the same reduced 
entropy in the thermodynamic limit n 



oo: 



,. lnW 2 , In Wi 
lim = hm 

n — >oo n n — >oo 



<^> lim — In 

n — >oo ri 



dG (I) 



01 



n 
SIq 



(500 



(9) 



which imposes the restriction that the Jacobian of the 
reparametrization can not grow exponentially with n. I 
shall assume for practical purposes the validity of the 
approximation S = Si = InWi ~ S2 = when the 

interest system is large enough. 
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III. OPEN SYSTEM 

Let us now consider an experimental setup where the 
interest system is put in contact with certain external ap- 
paratus prj which controls its macroscopic behavior in a 
way that the expectation values of the physical quantities 

= 9 (^fj are kept fixed: 



6' 



(10) 



The equilibrium distribution function <2i of the interest 
system under these conditions can be derived from the 
standard procedure of maximization of the entropy: 



S [<£>] = — Sp {uj IntD} , 



(11) 



which leads directly to the following generalized 
Boltzmann-Gibbs distribution: 



exp 



{- Vk e k (i)}. (12) 



The Lagrange multipliers rj = {rjk} are the canonical pa- 
rameters of the external apparatus p^, which will be re- 
ferred hereafter as a generalized thermostat or generalized 
reservoir. I shall show now that the above distribution 
function provides a suitable generalization of the well- 
known Boltzmann-Gibbs distributions: 



1 



Z{(3) 



exp 



(13) 



which is just a special case of the generalized distribution 
(JT2J) with 6 m = I. The convention A fe ■ tp k = A ■ ip will 
be hereafter adopted in order to simplify our notation. 
Since the operator: 



exp 



-77 -Of/ = / exp [—rj ■ (7)] 5 



I - I 



dl, 



(14) 

the canonical partition function Z (77) can expressed as 
follows: 



Z(ji) = J exp [-77 • 6 (I)] Cl (I) dl, 



(15) 



where Cl (I) = Sp |<5 I — I | is the microcanonical par- 
tition function in the representation TZi of the space 3. 

The integral (fTS"]) when (7) = 7 is just the Laplace 
transformation : 



Z(f3) = J exp [-/3 ■ I] Cl (7) dl, 



(16) 



which supports in the thermodynamic limit n —> 00 the 
validity of the well-known Legendre transformation: 



P((3)=M{(3-I-S(I)}, 



(17) 



between the Planck thermodynamic potential P ((3) = 
— In Z ((3) of the Boltzmann-Gibbs distribution (TTB"]) and 
the coarsed grained entropy S (I) = In [Cl (I) <57 ] of the 
microcanonical ensemble ([3]). The identity (j 1 7|) repre- 
sents that a macrostate of the Boltzmann-Gibbs ensem- 
ble (|13|) with canonical parameter f3 is equivalent in the 
thermodynamic limit 77, —> 00 to a macrostate of the mi- 
crocanonical ensemble Q whose set I* represents the 
point of global minimum of the functional: 



P{I;(3)=f3-I-S{I), 



(18) 



within the representation TZi of the abstract space 
which satisfies the conditions: 



d 2 S(P) 



(19) 



where the entropy Hessian Kij must be a negative defined 
matrix, that is, the entropy must be locally concave in 
the point I* E IZi. Such exigency implies that all those 
regions of the Euclidean subset Hi where the entropy 
S is not locally concave will never be equivalent to a 
macrostate of the Boltzmann-Gibbs ensemble (fl~3|) . This 
situation is usually referred as an ensemble inequivalence. 

Let Bj be the Euclidean n-dimensional subset com- 
posed by all admissible values of the canonical parame- 
ters [3 of the external apparatus leading to a thermody- 
namic equilibrium of the interest system characterized by 
the Boltzmann-Gibbs distribution (fT3"|) . According to the 
condition (|19[) , such thermodynamic equilibrium is char- 
acterized by the identification /3 = (3 of the canonical 
parameters [3 of the external apparatus with the com- 
ponents (3 of the gradient of the microcanonical entropy, 
the map ipj: 



0i = 



dS{I) 
dP 



(20) 



in the representation IZi of the space 9. When the en- 
tropy is locally concave in certain subset 7r Q C IZi , the 
map tpi establishes a bijective correspondence between 
the subset n a and certain subset A a C Bi. During the 
ensemble inequivalence the bijective character of the map 
tp is lost. The set (3 G Bi of canonical parameters of the 
external apparatus leads to a macrostate with ensemble 
inequivalence when there exist at least two distinct points 
I\ and I2 G Hi, where (3 = ipi = ipx {h)- Thus, the 
ensemble inequivalence can be seen as a thermodynamic 
equilibrium where exist a competition among all those 
microcanonical states Ik 6 Hi where ipi (Ik) = (3. 

Let us return to the case of the generalized Boltzmann- 
Gibbs ensemble (|12p . Taking into account the 
reparametrization invariance of the measure d/i = 
CI (7) dl = CI (0) d<d, the integral (fT5"l) can be rephrased 
as follows 



Z(r}) = J exp [-77 • 0] CI (0) dB, 



(21) 



which allows a natural extension of the conventional Sta- 
tistical Mechanics and Thermodynamics with a simple 
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reparametrization change I — > and (3 — > rj: The Leg- 
endre transformation between the thermodynamic poten- 
tials P (77) = - InZ (»j) and S (6) = In [Q (6) $0 O ]: 



inflri-e 



5(e)}, 



and the conditions of the ensemble equivalence: 

as(e*) d 2 s{&*) 



d® 



ft A: 



h kn 



(22) 



(23) 



where the entropy Hessian Kk m in the representation 1Z& 
must be a negative definite matrix, that is, the entropy 
S (0) must be locally concave in the point 9* € TZq. The 
above results deserve some remarks. 

Remark 1 While a reparametrization change does not 
alter the microcanonical description of an isolated sys- 
tem, for an open system such transformation implies 
a substitution of the external apparatus pp — > p r) 
which controls its thermodynamic equilibrium, and conse- 
quently, generalized Boltzmann-Gibbs ensembles with dif- 
ferent representation TZe 1 and !Ze 2 provide in general 
different thermodynamical descriptions. 

Remark 2 The validity of the Legendre transformation 
i22jl implies the local equivalence of the macrostates of 
different generalized Boltzmann-Gibbs ensembles in the 
thermodynamic limit n — > 00 whenever they are equiva- 
lent to the same macrostate of the microcanonical ensem- 
ble. Such equivalence implies a relation among the canon- 
ical parameters of the generalized thermostat and an or- 
dinary thermostat: Since S (I) = S (O) when n — > 00, 
the conditions U9\) and lead the following transfor- 
mation rule during a reparametrization change 0^1: 



ft = dS/dP 

rtk = ds/oe k 



a de k 
m = -op 1 ** 



while the entropy Hessian obeys: 

QQk QQW, Q2Qk 



dP dP 



dPdP 



(24) 



(25) 



Since IZi and TZq are just two admissible representa- 
tions of the abstract space 3, the above transformation 
rules are also applicable to any two admissible represen- 
tations of 3. Equation (124]) is just the transformation 
rule of the components of a covariant vector within a dif- 
ferential geometry. However, the entropy Hessian does 
not represent a second rank covariant tensor, since the 
transformation rule should be given by: 



QQk qqt, 

~dP~dP 



■Tkr, 



(26) 



instead of the rule 
third remark. 



This result leads to an important 



Remark 3 Since the entropy Hessian is not a second 
rank covariant tensor, the concavity properties of the en- 
tropy may change during the reparametrization changes. 
Consequently, the regions of ensemble inequivalence be- 
tween the Boltzmann-Gibbs ensemble and the micro- 
canonical ensemble in the representation IZi does not co- 
incide to the regions of ensemble inequivalence in the rep- 
resentation TZq . 

Let us consider a trivial example taken from our previ- 
ous work [T3|. Let s be a positive real map defined on a 
seminfinite Euclidean line C, s : C — > R + , which is given 
by the concave function s (x) = y/x in the representation 
1Z X of C (where x > 0). Let <p be a reparametrization 
change ip : 1Z X — > 1Z V given by y = <p(x) = i<. The 
map s in the new representation lZ y of the seminfinite 
Euclidean line C is now given by the function s (y) = y 2 
(with y > 0), which is clearly a convex function. 

The third remark implies a revision of our conceptions 
about the classification of the phase transitions based on 
the concavity of the entropy and the generalization of 
some Monte Carlo methods with distribution functions 
inspired on the Statistical Mechanics. 

The first-order phase transitions are closely related 
with the phenomenon of ensemble inequivalence. Since 
the ensemble inequivalence depends on the nature of 
the external apparatus controlling the thermodynamic 
equilibrium of the interest system, the existence of this 
anomaly reflects the inability of the external apparatus to 
control all microcanonically admissible equilibrium states 
of the interest system, and consequently, this kind of phe- 
nomenon can not be relevant within the microcanonical 
description. The interested reader can see a long discus- 
sion about this subject in the ref . (l7| . 

The phenomenon of ensemble inequivalence also im- 
plies an important the lost of information about the ther- 
modynamical properties of a given system during the 
phase coexistence in regard to the thermodynamical in- 
formation which can be derived from the microcanonical 
description. Such lost of information also leads to the 
failure of ordinary Monte Carlo methods based on the 
Boltzmann-Gibbs ensemble in the neighborhood of the 
first-order phase transitions [2(| ■ 

The ensemble inequivalence can be successfully 
avoided by using an appropriate representation TZq 
within the generalized Boltzmann-Gibbs distribution 
(fT2"|) . Therefore, the simple consideration of the 
reparametrization changes allows a suitable extension of 
some Monte Carlo methods. This question was discussed 
in the ref. [l8| for the particular case of the Gibbs canon- 
ical ensemble. I shall consider below a formal extension 
of the methodology developed in that reference for en- 
hancing the well-known Metropolis importance sample 
algorithm [2l| by focussing only the most important fea- 
tures. The reader could see the ref. [3 for more details. 

The probability of acceptance of a Metropolis move 
based on the generalized Boltzmann-Gibbs distribution 
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(Tp2f can be given by: 
P (I\I + AI) 



min {1, exp [—77 ■ A6]} , (27) 



where A8 fc = <d k {I + AI) 
enough AI « I => AQ k 
consequently: 



(I) . When n is large 
{de k (I) /dP) AP, and 



p(I\I + AI) = min jl, exp -/3 ■ AI X, 



(28) 



where the effective canonical parameter 0$ = fa (I; rf) is 
given by: 



dQ k (I) 



(29) 



Reader can notice that the above definition is very 
similar to the transformation rule l|24p. The difference 
resides now in the fact that the set I does not take the 
equilibrium value /*, but it changes during the Monte 
Carlo dynamics. Thus, the present generalized Metropo- 
lis algorithm looks-like an ordinary Metropolis algorithm 
where the canonical parameters (3 = [3 (I; rf) of the ex- 
ternal reservoir exhibits correlated fluctuations with the 
fluctuations of the physical quantities I characterizing 
the interest system around the corresponding equilibrium 
values (3 and I* connected by the ensemble equivalence 
within the generalized Boltzmann-Gibbs description (112|) 
and the transformation rule (|24|) . As already shown in 
ref. [l8| , this feature of the present Metropolis algorithm 
allows the Monte Carlo dynamics to explore the anoma- 
lous regions of the subset TZi where the ensemble in- 
equivalence within the Boltzmann-Gibbs description (fTS"]) 
takes place. The reader can see an example of application 
of this method in the section [V] The above observations 
lead directly to a fourth remark which clarifies us the 
nature of the generalized thermostat p n . 

Remark 4 When the interest system is large enough, 
the generalized thermostat or reservoir p v leading to the 
generalized Boltzmann-Gibbs distribution is almost 

equivalent to an ordinary thermostat pp whose canonical 

parameters (3 fluctuates with the fluctuations of the phys- 
ical quantities I of the interest system according to the 
rule fH. 



The fourth remark is particularly interesting. Conven- 
tional Thermodynamics usually deals with equilibrium 
situations where the fundamental physical quantities I 
are kept fixed when the interest system is isolated (a 
microcanonical description) or the physical quantities I 
fluctuate and corresponding canonical parameters [3 are 
kept fixed by the thermostat ("within canonical descrip- 
tion"). The generalized Boltzmann-Gibbs ensemble (fl2)) 
is just a natural framework for consider all those equi- 
librium situations where both the fundamental physical 
quantities I of the interest system and the corresponding 
effective canonical parameters (3 of the external thermo- 
stat fluctuate as a consequence of their mutual interac- 
tion. Thus, the ordinary condition for the thermody- 
namic equilibrium between the canonical parameters of 



the interest system (3 and the external apparatus (3 only 
takes place for their corresponding average values: 











(30) 



The last observation leads directly to the main objective 
of the present study: the implementation of a general- 
ized fluctuation theory based on the reparametrization 
invariance of the microcanonical ensemble. 



IV. GENERALIZED FLUCTUATION THEORY 
A. Fundamental results 

Let 51 = I — I* and 5/3 = (3 — j3 be dispersions of the 
fundamental physical quantities I of the interest system 
and the effective canonical parameters (3 of the general- 
ized thermostat respectively. Our aim is to obtain the 
correlations (5P6P), (5/3i5P) and (5f3iSf3j) within the 
generalized Boltzmann-Gibbs description (TT2"|) and iden- 
tify their mutual relationships. It is very important to 
remark that such correlations depends crucially on the 
interaction between the interest system and the gener- 
alized thermostat, and consequently, the values of these 
quantities are modified with the change of the external 
control apparatus. 

Let us began from the correlation (SPSP): 

(SP5P) = — ^- [ 5P5P exp \-r) • 6 (1)1 fi (I) dl, 

Zyn) J 

(31) 

where the Gaussian approximations leads directly to the 
result: 



(5P5P) = (njf 

where the matrix IK is given by: 



(I) = 9 2 9 (J*) _ d 2 S(P) 
ij gjigjj gjiQjj ■ 



(32) 



(33) 



The above result generalizes the classical relation among 
the correlations (5P5P) and the entropy Hessian 
(TTl]), (5P5P) = - Ky 1 derived fr om the Boltzmann- 
Gibbs distribution (TL3")) . Taking into consideration the 
transformation rule (|25p . the matrix H^y can be rewrit- 
ten as follows: 



n 



(-0 _ 



dP OP 



(34) 



where n^m is the entropy Hessian (|23[) in the representa- 
tion TZq of the abstract space 3. 

The expectation values 6* = (6 fc ) and (50 k 5& m ) 
with (50 = G — 0* within the generalized Boltzmann- 
Gibbs description (TT2")) can be obtained from the Planck 
thermodynamic potential P (rf) in the representation TZq 
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in complete analogy with the standard Statistical Me- 
chanics: 



Finally, the approximation (|40|) and the identity (|4T1) 
allows us to express the correlations (SfiiSftj) as follows: 



(-)- .(a*- 



(se k se m ) = - 



d 2 P (??) 
dr] k dr] m 



(35) 



The Gaussian estimation of (<5O fe <50 m ) is given now by 
the inverse of the entropy Hessian 



(36) 



(^e m ) = -^ = (n^) \ 

Taking into account the result (|34|) : 



n 



(f) _ae*fle"- n( e, 



si* a# 



km ' 



(37) 



and therefore, the correlations (SPSP) and ^<50 fc (5O' m ) 
are related by a transformation rule of second rank con- 
travariant tensors: 

where the correlations (SPSP) can be expresses in terms 
of the Planck thermodynamic potential P (if) as follows: 



(SPSP) 



dP dP d 2 P(ri) 
~gQk d Qm dr) k dr) m 



(39) 



The correlations (SPiSP) and (Sj3iS(3j) can be derived 
easily from the correlations (SPSP). Since SI << I* 
when the interest system is large enough, the dispersion 
8(3 of the effective canonical parameters = j3 (I; rf) de- 
fined in the equation (|29p can be estimated in terms of 
the dispersions SI as follows: 



d 2 & k (I*) 



(40) 



and therefore: 



Since the definition 



of the matrix IT- J ' allows us to 



(41) 



express the matrix f ^ derived from the second deriva- 
tives of the map as follows: 

fl»8*(J*) n (7) , d*S(I*) 
the correlations (Sf3iSP\ can be given by: 

(JAM*) = Sf + (SISP) , (43) 

where the delta Kronecker Sf appears as consequence of 

the presence of the terms II^p (SI n SP) = Sf which is 
straightforwardly followed from the equation 



(SPiSfy) = F in F jm (SI n SI m ) = F jm (SfrSI™) , (44) 
and considering the relation: 

F jm E(^if)IlH (45) 
I arrive to the following result: 

{8^) = (Sftsn n« (Sfosn . (46) 

The relations (|43|) and (|46|) are the fundamental re- 
sults of the present subsection. Introducing the nomen- 
clature for the correlations tensors G y = (SPSP), M\ = 
(SPiSP) and Tij = (SfliSflj), the fundamentals relations 
can be rewritten as follows: 



Mi=Sl + K im G m ^ 



— J-J-nrr, 



(47) 
(48) 



where G nm = (G l ) 

I can substitute (|47|) in (|48l in order to express the 
correlation matrix Mf in terms of the entropy Hessian 
Ha and the correlation matrix G lJ : 



= M™ (G m + Ki„) , 

= (<5? + Kjm G mn ) (G m + «<„) , 

which can be rewritten as follows: 



(49) 



G — Sj + 2G + G Ki n G Kj m 
= (SI + G m n l0 ) («jf + G^ m ) , 



(50) 



Therefore, the identity (|48|) can be expressed as follows 



(51) 



The reader may notice that no one of the above iden- 
tities makes any explicit reference to the representation 
IZq used in the generalized Boltzmann-Gibbs ensemble 
Q12[) . and consequently, such thermodynamical identities 
possesses a general validity for the case of an external ap- 
paratus controlling the average values of the fundamental 
quantities I of the interest system throughout the iden- 
tification of the average values of the effective canonical 
parameters (3 of the external apparatus with the corre- 
sponding average values of the parameters (3 — dS/dl 
derived from the gradient of the microcanonical entropy 
S of the interest system. 

Since the thermodynamical variables I and [3 are just 
an admissible representation from the reparametrization 



6 



invariance viewpoint, the reparametrization change I - 
O and (3 — > T) leads also to the validity of the relations: 



d 2 S(Q) 



k de k oe 
(se r se m ) = (s Vk SQ n ) (5 Vn se r ) , 



(se n se m ), (52) 

(53) 



which means that the fundamental results also exhibit 
reparametrization invariance. 

B. Some anomalous behaviors revised 

As elsewhere shown [22[ , the entropy Hessian provides 
important information about the fluctuations of thermo- 
dynamic quantities during the thermodynamic equilib- 
rium of an open system. Such connection is usually de- 
rived from the Boltzmann-Gibbs distribution (|13|) which 
leads directly within the Gaussian approximation to the 
following result: 



d 2 s 

dPdP 



-K im (6I m SP)=6i. (54) 



The correlation matrix G 1 -? = (6P6P) is always nonneg- 
ative, which is easily verified from the following reason- 
ings: 

((^if) = &6 (5P5P) = > 0, (55) 

where the equality only takes place when £ = 0. Conse- 
quently the applicability of the relation ([54")) is limited to 
those regions of the subset IZi where the entropy Hessian 
be a negative definite matrix, that is, the regions in which 
is ensured the equivalence between the Boltzmann-Gibbs 
ensemble (|13[) and the microcanonical ensemble. Thus, 
the extrapolation of the formula (|5"4"|) towards regions of 
ensemble inequivalence has no sense within the conven- 
tional Thermodynamics. Let us consider a very simple 
example. 

Let be a system whose thermodynamical state is de- 
termined only by the total energy E. The corresponding 
canonical parameter is the inverse temperature (3 of the 
Gibbs thermostat. The relation 154[) is given now by: 



(56) 



\ / ~ 1' 



which can be conveniently rewritten by using the mi- 
crocanonical inverse temperature (3(E) = 1/T(E) = 
dS (E) /dE as follows: 



(dP{E)Y 1 
V dE 



,dE 
dT 



= T 2 C, 



(57) 



where C is the heat capacity. Since the average square 
dispersion of the total energy energy (SE 2 ^ > 0, the 



heat capacity should be nonnegative C > 0. As already 
mentioned, such relation is only applicable during the 
ensemble inequivalence. 

There is nothing wrong in assuming that all those en- 
ergetic regions with a convex entropy d 2 S (E) /dE 2 > 
are just thermodynamical states with a negative heat ca- 
pacity C < 0. Such anomalous behaviors are well-known 
since the famous Lyndel-Bell work [16j | and them have 
been also observed recently in experiments involving nu- 
clear or molecular clusters [H, HH [HJ • However, the 
relation between the heat capacity C with the average 
energy fluctuations (|57[) has to be generalized in order 
to attribute some reasonable physical meaning to these 
anomalous thermodynamical states in terms of a suit- 
able fluctuation theory. Such generalization is provided 
precisely by the formula: 



(5/36E) 



1 



d 2 S(E) 
dE 2 



(5E 2 ) 



(58) 



which is a particular expression of the fundamental result 
(|43p to the unidimensional case (only one control param- 
eter). While the inverse temperature can be kept fixed 
5(3 = in order to control a microcanonical state with 
a positive heat capacity C > 0, such control parameter 
can not be kept fixed when we deal with microcanonical 
states with non positive heat capacities C < 0. There 
(8(35E) > 1, which means that the fluctuations (3 and / 
should be correlated in order to ensure the control of such 
anomalous thermodynamic states, that is, in order to en- 
sure the ensemble equivalence. As already commented in 
the subsection above, this feature is the key for the suc- 
cess of the Metropolis algorithm based on the generalized 
Boltzmann-Gibbs ensemble (Tl"2"|) . 

The heat capacity C is a very simple example of a re- 
sponse function which characterizes how sensible is the 
system energy E (controlled quantity) during a small 
change of the temperature T of the thermostat (control 
parameter) . Another example of response function is the 
magnetic susceptibility \b = dM/dB, which character- 
izes how sensible is the system magnetization M (con- 
trolled quantity) during a small change of the external 
magnetic field B (control parameter). Thus, the well- 
known thermodynamical identities: 



C 



dE 
dT 



(3 2 (5E 2 ) , X b 



dM 
~dB 



(3(5M 2 ), (59) 



can be considered as response-fluctuation theorems. As 
already pointed out, such relations are non applicable 
during the ensemble equivalence where C or xb can 
admits negative values, and therefore, such thermody- 
namic theorems should find a natural extension within 
the present framework. 



C. Microcanonical thermodynamic formalism 

The geometric framework developed in the present 
work constitutes a natural generalization of the thermo- 
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dynamic formalism of the conventional Thermodynam- 
ics. A key for the success is to perform a thermodynam- 
ical description whose essence be so close to the macro- 
scopic picture provided by the microcanonical ensemble, 
which possesses a hierarchical supremacy in regard to any 
generalized Boltzmann-Gibbs description (fT2|) . 

The total differentiation of the microcanonical entropy: 



dS = fad! 1 with ft 



dS 



(60) 



provides the expectation values of the effective canonical 
parameters (3 of an external generalized thermostat and 
the physical quantities / in the representation IZi of the 
abstract space 9. The symmetry of the entropy Hessian 
Kij leads to what could be considered as the Maxwell 
identities within the microcanonical description: 



d 2 S 
dPdP 



ah 

dl 



The reparametrization invariance imposes some re- 
strictions to the admissible parameters used in the gener- 
alized Boltzmann-Gibbs canonical description within the 
present geometric framework. As already shown in the 
equation (|24|) . the canonical parameters ft are just the 
components of a covariant vector. Such status can not be 
attributed to the ordinary temperature T or the external 
magnetic field B. For example, the Gibbs canonical en- 
semble written for a ferromagnetic system with Hamilto- 
nian H under an external magnetic field B directed along 
the z-axis: 



1 



Z(P;B Z ) 



exp 



H - BM Z 



(62) 



can be considered as a Boltzmann-Gibbs distribution 
with canonical parameters (3 and A = —f3B, whose mi- 
crocanonical description is given by: 



1 



tt(U,M z ) 



5[U-H)5[M. 



M, 



(63) 



While the parametrization (ft A) is just a covariant vector 
within a differential geometry, such character is not ex- 
hibit by the the pair (T, B), and therefore, this last is not 
a suitable parametrization of the generalized canonical 
description from the reparametrization invariance point 
of view. The above observation inspires a reasonable and 
necessary redefinition of the concept of response functions 
in order to preserve the geometric covariance of the ther- 
modynamic formalism developed in this work. 

The element of the generalized response matrix \ lJ '■ 



d(r) 



(64) 



characterizes the rate of change of the average of the i- 
th physical quantity (P^ under a small variation of the 



j-th average effective control parameter ft 



ft of 



the generalized thermostat. Thus, x y is a measure of 
the sensibility of the interest system under the control of 
the external apparatus. Strictly speaking, the response 
matrix x lJ (|64p is not a second-rank contravariant tensor. 
Such character is straightforwardly followed from the fact 
that the microcanonical counterpart Xij °f the inverse 
response matrix Xij- 



Xij 



_9ft 
d(P 



is just the negative entropy Hessian Xij = ~ K ij UM 
which obviously is not a second-rank covariant tensor. 
The connection between the generalized response ma- 
trix x u with the entropy Hessian K.y leads directly to 
rephrased the microcanonical Maxwell identities (|6"Tj) as 
the symmetric character of the generalized response ma- 
trix: 



(65) 



The ensemble equivalence within the generalized 
Boltzmann-Gibbs ensemble (|12l) in the representation TZi 
demands the negative definition of the entropy Hessian 

, which is equivalent to demand the positive definition 
of the generalized response matrix x*- 7 (|6"4"|) . Thus, regions 
of ensemble equivalence are just regions with anomalous 
behavior in the response functions x lJ ■ 

The correlations matrixes G y = [8P8P), M'( = 
(SftiSP) and Ty = (JftiJft) are related by a set of ther- 
modynamical identities which could be referred as the 
generalized fluctuation relations: 



(S&5P) 
(SP6I n ) (60 n 5l3j) 



8\ + Kij(8I m 8P), (66) 
(SftSI") (S/3JP) , (67) 



within the representation IZi. The connection between 
the entropy Hessian Kij with the generalized response 
matrix x tJ allows us to rephrase the first-generalized re- 
lation (RHH) as follows: 



X ii =x im (Sf3 m 8P) + (5PSP) 



(68) 



which is the generalized response-fluctuation theorem. 

A fundamental difference of the above expressions in 
regard to the conventional result ([54")) is the presence 
of a non vanishing correlation matrix M? — (8{3i8P\. 
Moreover, the identity (|6T[) has no counterpart within 
the conventional Thermodynamics. Such correlations are 
not present within the ordinary Boltzmann-Gibbs distri- 
bution (|13[) since the canonical parameters of the ther- 
mostat are keep fixed, so that 8(3 = (S/3i8P) = 
(8f3 n 8f3j) = 0. Such limitation does not allow to the 
Boltzmann-Gibbs description to access to all those re- 
gions with a anomalous behavior of the response func- 
tions, i.e. regions with a negative heat capacity. As al- 
ready shown by us in the refs.JU EH, the only way to 
access to such anomalous regions with ensemble inequiva- 
lence within a Boltzmann-Gibbs-like description is by us- 
ing a generalized thermostat with a fluctuating effective 
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canonical parameters (3 = r\ • dQ/dl in order to control 
the thermodynamic equilibrium of the interest system by 
keeping fixed the averages of the ordinary equilibrium 

conditions (fiij — (Pi) ■ 

The generalized fluctuation relations (|66[) admit an- 
other interpretation. Let us consider a representation 
TZj where the entropy Hessian be a diagonal matrix in 
a given point J, k%$ — K%5i$, which can be obtained 
from the representation IZi by using an appropriate 
orthogonal rotation. Whenever an Hessian eigenvalue 
Ki > (when the interest point belongs to a regions 
with ensemble inequivalence) the corresponding correla- 
tion (SfliSJ 1 } = 1 + Ki ^(<5J*) 2 ^, and consequently: 

(SPiSr) > 1. (69) 

Such inequality looks-like an uncertainly relation which 
talks about that the physical quantity J 1 and the corre- 
sponding canonical parameter /?$ of the generalized ther- 
mostat can not be kept fixed when the corresponding 
eigenvalue m of the entropy Hessian be nonnegative, that 
is, when the thermodynamical state of the interest system 
belongs to a region of ensemble inequivalence. Thus, the 
quantities /3j and J J can be considered within the present 
framework as complementary thermodynamic quantities. 
This interpretation allows us to refer the generalized fluc- 
tuation identities (|66p as uncertainly thermodynamic re- 
lations pTj . 

It is necessary before end this subsection to devote 
some comments about the phase transitions. The phase 
transitions within the conventional Thermodynamics are 
mathematically identified by the non analyticities of the 
thermodynamic potential which is relevant in a given ap- 
plication. Within the microcanonical description of an 
isolated system the only relevant thermodynamic poten- 
tial is the microcanonical entropy S. 

The microcanonical entropy S is always a continuous 
function, and it is analytical whenever the number of de- 
grees of freedom of the interest system n be finite. Thus, 
as in the conventional Thermodynamics, the non analyt- 
icities of the entropy per particle s = S/n can appear only 
in the thermodynamic limit n — > oo. Since the micro- 
canonical entropy is a scalar function from the viewpoint 
of the reparametrization invarian.ee, its nonanalyticities 
will appears in any representation TZq of the abstract 
space Sy, and therefore, these mathematical anomalies are 
reparametrization invariant. Thus, the lost of analyticity 
of the microcanonical entropy is identified in the present 
geometric framework as the mathematical signature of 
the microcanonical phase transitions. It is said micro- 
canonical phase transitions in order to distinguish them 
from other macroscopic anomalies which are related with 
the occurrence of the phenomenon of ensemble inequiv- 
alence in an open system. Generally speaking, there ex- 
ist two generic types of microcanonical phase transitions: 
the microcanonical discontinuous phase transitions (asso- 
ciated to a discontinuity in some of the first derivatives 
of the microcanonical entropy S) and the microcanonical 



continuous phase transitions (associated to lost of ana- 
lyticity of the microcanonical entropy S with continuous 
first derivatives). 

The dynamical origin of the microcanonical ensemble 
leads to identify any microcanonical phase transition as 
the occurrence of some kind of dynamical anomaly in 
the microscopic picture of the system. The evidences an- 
alyzed in the previous work [ItJ suggest strongly that the 
microcanonical phase transitions are directly related with 
the occurrence of the ergodicity breaking phenomenon, 
which takes place when the time averages and the ensem- 
ble averages of certain macroscopic observables can not 
be identified due to the microscopic dynamics is effec- 
tively trapped in different subsets of the configurational 
or phase space during the imposition of the thermody- 
namic limit n — > oo [23| ■ 

As elsewhere shown, the convexity of the microcanon- 
ical entropy leads to the existence of a lost of analyt- 
icity of the Planck thermodynamical potential P ((3) in 
the thermodynamic limit via the Legendre transforma- 
tion (|17p . Such discontinuity in some of the first deriva- 
tives of P (f3) is the signature of the first-order phase 
transitions in the conventional Thermodynamics. It is 
necessary to remark at this point that this kind of anoma- 
lies can not be microcanonically relevant since they have 
only sense for the open systems instead of the isolated 
ones [l7| . The first-order phase transitions are untimely 
related with the phenomenon of ensemble inequivalence 
appearing as a consequence of the non negative definition 
of the entropy Hessian k^- in a given representation TZq 
of the abstract space Therefore, the existence of the 
ensemble inequivalence depends crucially on the nature 
of the external apparatus controlling the thermodynamic 
equilibrium of the interest system. Thus, the first-order 
phase transitions can be considered as avoidable ther- 
modynamical anomalies of the open systems. This last 
feature can be used for enhancing the available Monte 
Carlo methods based on the Statistical Mechanics in or- 
der to avoid the ensemble inequivalence by considering 
the generalized Boltzmann-Gibbs distributions (|12p. 

I recommend to the interested reader to read the long 
discussion exposed in the ref. [ijj for a better understand- 
ing about all these important questions. I shall consider 
in the section[V]ari example where the existence of a first- 
order phase transition in the canonical description hides 
the occurrence of two microcanonical phase transitions 
in the thermodynamical description of the q = 10 states 
Potts model system. 

D. Conjugate variables and the Riemannian inter- 
pretation of the fluctuations 

The introduction of the dispersion Scot : 

5u>i = 5f3 l - K m 5I n , (70) 

allows to rewrite the fundamental relation as follows: 

(6^51^ = 6?. (71) 
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The correlations (SuiiSuij) can be derived as follows: 

= (SfiiSiOj) - K in (8I n 5u)j) , 



The Gaussian approximation allows us to estimate the 
canonical partition function as follows: 



Z (77) = cxp [-P*] ^det(27rGy )<5/ - 1 , 



(81) 



(72) 



where the substitution of equation (|47|) and (|49j) leads 
finally to the following result: 



which allows us to estimate finally the generalized 
Boltzmann-Gibbs distribution in terms of the variable 
I by: 



(SuiSujj) = Gij. 



(73) 



duj GC {1} = 



1 



v /det(27rG i J) 



cxi. <j —GijSPSF } ill. 



It is very easy to notice that the quantity u>i is just the 
difference between the effective canonical parameter /3 = 
rj ■ <99 I dl of the external apparatus p v and the canonical 
parameter of the system (3 = dS/dl, Wj = ,9,-/3,, that 
is: 

d& k dS d , n , 
u Jl = r lk —-— = —P(Q-r l ), (74) 

where P (9; 77) is the functional: 

P(e;r?) =r)-Q-S. (75) 

The quantities to = {u>i} and I = {/*} can be classified 
as conjugated variables since they obey the relations: 



The transformation from / towards the conjugated vari- 
ables uj allows us to rewrite the above results as follows: 



duJcc W} = duj J S {lu — lu (I)} duJac {1} ■ 



(83) 



which leads immediately within the Gaussian approxi- 
mation to the result: 



du>GC {^} 



y/det (2irG t 



: exp 



1 



G^uj.ujj \ duj. (84) 



The correlation tensor Gij or G 13 arises as a conse- 
quence of the interaction between the generalized ther- 
mostat (Fij) and the interest system (k^) within the gen- 



(5PSP) = G ij , (SuiSr) = (SutSuj) = G i:j . (76) eralized Boltzmann-Gibbs description {G i3 = Fj 



In fact, the approximation (JJDJ) leads directly to the fol- 
lowing transformation rule between the fluctuations 5u>'s 
and Si's: 



Such tensor provides in this context a Riemannian inter- 
pretation of the fluctuations, where the norm: 



Ss 2 = GuSrSP, 



(85) 



8 = F ij SP ^Suii = Gij5F, 



(77) 



which allows us to obtain the second and third identities 
from the definition of the correlation matrix G 13 in the 
equations ([76]) . 

The generalized Boltzmann-Gibbs distribution (TT2|) 
can be rephrased as follows: 



dUGC {1} = ^TT CXP [-V ■ e (I) + S (I)} 



is directly related to the probability of occurrence of a 
given fluctuation 51 = I — P around the thermodynamic 
equilibrium of the system in accordance with the equa- 
tion (18211. I refer to a Riemannian metric since the cor- 



relation matrix behaves as a second-rank tensor un- 
der a given reparametrization v : IZi — ► IZj within the 
same representation TZq of the generalized Boltzmann- 
Gibbs ensemble. Taking into account J = v (7 ) =>■ 
9 (I) — 9 [v^ 1 (J)] , the expression of the correlation 
matrix in the representation IZj : 



1 



exp[-P(Q;ri)}dI. (78) 



G a @ 



The power expansion of the functional P(Q;rj) in the 
neighborhood of its minimum is given by: 

P (9; V )=P (9*; 77) + dP ^*J^ §p 



d 2 P(e*;n) 

dJ a dJP 
dP dP 



dP dP 
d.J a dJP ij 



d 2 P 



dJ a dJP 



(86) 



ld 2 P{e*;n) 
2 dPdP 



dP 
5PSP 



= P* + cu*SP + -GijSPSF 



(79) 



where 



dP(Q*; V ) _ _ d 2 P(Q*;n) 
- — - U, G i3 - aTiaTi ■ (8U) 



dP 



dPdP 



where the vanishing of the second terms takes place as a 
consequence of the stationary conditions (|80p . This re- 
sult generalizes the other Riemannian interpretations of 
the Thermodynamics obtained in the past [13]. While 
the covariant second-rank tensor Gij characterizes the 
probabilities of the fluctuations of the physical quanti- 
ties I = of the interest system, its corresponding 
contravariant tensor G %3 characterizes the probabilities 
of the fluctuations of the conjugated variables lu = {uJi} 
in the neighborhood of its expectation lu* = (lu) = 

0)— (/3) = in accordance with the equation l[8"3]). 
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which provides us information about the unusual thermo- 
dynamic equilibrium between the generalized thermostat 
and the interest system. Thus, the conjugated charac- 
ter of the variables u> and I is closely related with the 
complementary character of the canonical parameters (3 
or (3 and their corresponding physical quantities of the 
interest system /. 



V. A FORMAL APPLICATION 

Let us reconsider again the example about a ferromag- 
netic system under the influence of an external magnetic 
field B directed along the z-axis. Let us to introduce a 
distinction between the internal Hamiltonian H (related 
with the system internal energy U) and the total Hamil- 
tonian Hb = H — BM Z (related with the total energy 
E = U — BM Z ). As already commented, the Boltzmann- 
Gibbs distribution of the system under these conditions: 



1 



LOG 



Z(fcB z ) 



exp 



-PH E 



(87) 



can be reinterpreted as a Boltzmann-Gibbs distribution 
(BGE) with canonical parameters (3 and A = —f3B: 



1 



Z{{3;X) 



exp 



-f3H - XM Z 



(88) 



The above observation allows to distinguish the usual 
microcanonical description (IME): 



1 



% = 



n{E;B) 



5 E- H f 



(89) 



from the detailed microcanonical description (DME): 



1 



LOm = 



rt(u,M z 



6 (U - H 



M z - M, 



(90) 



While the canonical ensembles (|5T1) and ([55)) are iden- 
tical, the microcanonical ensembles (|89|) and (l90|) are 
essentially inequivalent. The above interpretation al- 
lows us to consider the modulus of the external mag- 
netic field B or more exactly, the quantity A = —f3B 
as the canonical control parameter which is thermody- 
namically complementary to the magnetization of the 
system M z . Consequently, the macroscopic description 
provided by ensemble (|89p can be considered as an in- 
termediate description between the DME (|9T)|) and the 
BGE (|88|) . which leads to the following hierarchy among 
these ensembles: DME=^IME=»BGE. This hierarchy pre- 
supposes that the description DME contains more infor- 
mation than the IME description, and at the same time 
this last one contains more information than the BGE 
description. 

Let us consider the analysis of the above magnetic sys- 
tem in terms of the thermodynamic formalism. The total 
differential of the microcanonical entropy of the detailed 
ensemble (|90p S (U, M z ) is simply given by: 



where: 



dS = f3dU + XdM, (91) 
8(3 



OX 



dU 1 dM z \dUj M \dM z/l 



(92) 



The introduction of the total energy U = E+BM Z allows 
to rewrite the above expression as follows: 

dS = (3 (dU - BdM) = [3{dE + M z dB) , (93) 

from which are derived the relations: 

S {E; B) = S (U, M z ) , (94) 



(95) 

Thus, the consideration of the total energy E instead of 
the internal energy U corresponds to a reparametrization 
where (U,M Z ) (E,B). Reader can notice that such 
transformation does not correspond to the diffcomorphic 
transformation considered by the reparametrization in- 
variance of the microcanonical ensemble discussed in the 
present work. Finally, the Planck potential P (f3 1 A) of 
the BGE is derived from the DME microcanonical en- 
tropy S (U, M z ) throughout the Legendre transformation 
as follows: 



P {(3, X)=(3U + XM Z - S (U, M z ) , 



(96) 



which is equivalent to the one obtained from the IME 
microcanonical entropy S (E; B) by considering the rela- 
tions A = -(3B and E = U - BM Z : 

p ((3, A) = f3U — (3BM - S (U, M z ) 

= (3E - S (E; B) = P (/?; B) . (97) 

The microcanonical partition functions f2 (U, M z ) and 
f2 (E; B) are related as follows: 



Q(E;B) = J Q (U, M z ) 8 (U 



E - BM Z ) dUdM z , 



J fl(E + BM z ,M z )dM z , 



(98) 



which leads to the following maximization problem in the 
thermodynamic limit: 

S (E; B) = sup {S(E + BM Z ,M Z )}, (99) 

whose stationary conditions are given by: 
dS(U,M z )„ , dS(U,M z ) 



dU 



-B 



B' 



,d 2 s 



2B 



dM z 
d 2 S d 2 S 



0, (100) 



dU 2 dU3M z DM? 



+ TTITo < 0. (101) 
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The relation (|99[) clarifies the precise meaning of the re- 
lation (|94p . which talks us that the reparametrization 
(U, M) — » [E, B) in the thermodynamic limit is just a 
kind of projection where it could be involved an impor- 
tant lost of thermodynamic information. 

Taking into consideration that j3 — dS/dU and A = 
dS/dMz, the condition (jlOOj) is just the relation A = 
—(3B. On the other hand, the general solution of this 
stationary condition leads to what could be called as the 
microcanonical magnetization curve at a given magnetic 
field B, M — M (U;B). Generally speaking, such micro- 
canonical magnetization curve could be a multivaluated 
function of the internal energy, that is, there could be 
more than one admissible value of the magnetization M z 
at a given value of the internal energy U . However, only 
those points of the microcanonical magnetization curve 
satisfying the stability condition (|10ip will be observed 
within the IME description. Such condition is related 
with the negative definition of the entropy Hessian Kij, 
which can be rephrased by eliminating the magnetic field 
B as follows: 

f = X 2 kuu + A/3 (k um + kmu) + kmmP 2 < 0, (102) 



1.60 



0.2 0.4 0.6 0.8 1.0 1.2 



where the stationary condition (j 1 00!) was taken into ac- 
count. Thus, all those points where £ > can not be 
observed in the IME description, and consequently, such 
possibility involves a lost of thermodynamic information 
in regard to the one provided by DME description. Thus, 
some branches of the microcanonical magnetization curve 
can disappear within the IME description, leading in this 
way to the existence of discontinuities. Such discontinu- 
ities can be considered as the signature of microcanonical 
discontinuous phase transitions within the IME descrip- 
tion because of them represent discontinuities in the first 
derivative 



(103) 



of the IME microcanonical entropy S (E; B). The lost of 
information is more significant in the BGE description, 
which demands now the negative definition of the entropy 
Hessian: 



kuu < 0, kmmkjju - kumkmu > 0, 



(104) 



which is more restrictive than the condition (|102p . that 
is, the attainability of the condition (|104|) leads to the 
attainability of (|102p , but the converse is not true. 

The picture described above can be observed during 
the study of the thermodynamic properties of the q = 10 
states Potts model [24T ): 



H 



£ 

(ij) 



l-6„ 



(105) 



where the sum involves only neighbor-neighbor interac- 
tions on a square lattice L x L with periodic bound- 
ary conditions. The spin variable of the i-th site u\ = 



i — < — i — ' — r 
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FIG. 1: Comparative study between the Metropolis Monte 
Carlo based on the generalized canonical ensemble (GCMMC) 
and the Swendsen-Wang (SW) cluster algorithm during the 
thermodynamical description of the q = 10 states Potts 
model. The SW algorithm is unable to describe the mi- 
crocanonical states with a negative heat capacity during the 
first-order phase transitions since this method is based on the 
consideration of the Gibbs canonical ensemble. This aim is 
sucessfully performed by using the GCMMC algorithm, which 
is able to explore all those energetic regions of ensemble in- 
equivalence. Notice the bimodal character of the energy dis- 
tribution function within the Gibbs canonical description at 
f3 = 1.42, which is a feature of the first-order phase transitions 
illustrating the competition among two metastable states with 
different energies (coexistence phases). 



1,2, ...10 are rephrased as bidimensional vector — 
[cos (wtTi) , sin (w<7i)] with lu = 2n/q, which allows to in- 
troduce the microscopic magnetization as follows M = 

The study of this model system with L — 25 by us- 
ing the Gibbs canonical ensemble ([57]) is characterized 
by the existence of a discontinuous phase transition at 
(3 C ~ 1.42 with B = 0, which is recognized by the 
existence of a plateau in the caloric curve (3 versus e 
(e = E/N is the energy per particle with N = L 2 ) 
from e\ = 0.319 to £2 — 1.099 corresponding to a la- 
tent heat qn = £2 — £1 — 0.78. As already commented, 
the first-order phase transition within the Gibbs canon- 
ical ensemble is a consequence of an ensemble inequiva- 
lence between this description and the IME description 
which takes place as a consequence of the existence 
of thermodynamical states with a negative specific heat. 
These anomalous states manifest themselves as a back- 
bending in the caloric curve within the IME microcanon- 
ical ensemble 
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FIG. 2: Magnetic properties of the q — 10 states Potts model 
at zero magnetic field B within the IME microcanonical de- 
scription in which is shown the existence of two microcanon- 
ical phase transitions at s/f ~ 0.7 (ferro-ferro, discontinuous 
PT) and £f P — 0.8 (ferro-para, continuous PT). The dash-dot 
lines indicate the limit of the energetic region with a negative 
heat capacity hidden by the ensemble inequivalence. 



Reader can see all these behaviors in the FIG[TJ which 
shows a comparative study between two Monte Carlo 
simulations of this toy system at B = 0: the well- 
known Swedsen-Wang clusters algorithm (SW) based on 
the Gibbs canonical ensemble [2fj] and the Metropolis 
Monte Carlo algorithm (GCMMC) based on the gener- 
alized Boltzmann-Gibbs ensemble (IT2t explained in the 
section IIIII and the ref . [l8( . The S W algorithm is unable 
to describe all those states with a negative heat capacity, 
a task successfully performed by using the GCMMC al- 
gorithm. The key for the success of this last Monte Carlo 
method is precisely the using of a generalized Gibbs ther- 
mostat with a fluctuating inverse temperature (3. The 
caloric curve is obtained here by computing the average 

values (3 = ds(s;B) jde — (j3) and e = (e), while the 

curvature or entropy Hessian n = d 2 s (e; B) /de 2 can be 
derived from the unidimensional generalized fluctuation 
relation: 



(606e) = 1 + k (Se 2 ) 



by computing the correlations (5/35e) and (<5e 2 )- Reader 
can notice that no abrupt change of the energy appear 
by using the GCMMC algorithm, which illustrates us the 
avoidable character of the first-order phase transitions, 
and therefore, the irrelevance of such anomalies within a 
microcanonical description. 

The study of the magnetic properties of this system 



within the IME microcanonical ensemble (|89| at B = 
is shown in the FIGjH This figure reveals the existence 
of two microcanonical phase transitions which are hidden 
by the ensemble inequivalence within the Gibbs canon- 
ical description. It is easy to show that such anomalies 
are directly related with the lost of analyticity of the 
microcanonical entropy s (e; B) in the thermodynamic 
limit as well as the occurrence of ergodicity breaking 
in the microscopic picture of the system. The inter- 
ested reader can see more details about this study in the 
refs.JT?], EH- I n spite of the significant amount of ther- 
modynamical information provided by the IME descrip- 
tion (|89|) in regard to the BGE l[8?|). the discontinuous 
character of the magnetization curve is a clear indica- 
tion about the lost of information of the IME description 
in regard to the DME one (|90p appearing as a conse- 
quence of the unattainability of the condition (| 1 02|) for 
some macrostates. The Monte Carlo simulations shows 
the existence of several metastable states with differ- 
ent magnetizations in the energetic region (ea,£b) with 
ea — 0.7 and £b — 0.93, which indicates clearly the mul- 
tivaluate character of the microcanonical magnetization 
curve M = M (U; B = 0) in the DME description, which 
should exhibit an anomalous branch characterized by an 
increasing of the magnetization with the increasing of the 
internal energy dM/dU > 0. 

In complete analogy with the Gibbs canonical ensem- 
ble case, the existence of metastable states provokes an 
exponential divergence of the correlation times with the 
increasing of the system size during the Monte Carlo sim- 
ulations, a dynamical phenomenon referred in the litera- 
ture as supercritical slowing down. The existence of this 
phenomenon is untimely related with the lost of informa- 
tion associated to the ensemble inequivalence. To avoid 
this difficulty is necessary to perform a thermodynamic 
description within the DME ensemble. This aim could 
be carried out in an appropriated canonical-way within 
a generalized Boltzmann-Gibbs description (|12[) whose 
generalized thermostat exhibits two fluctuating effective 
control parameters (3 and A. This is just an unusual 
equilibrium situation from the conventional Thermody- 
namics viewpoint characterized by a fluctuating temper- 
ature T = 1/(3 and a fluctuating external magnetic field 
B = -\/[3. 

The dependence of the canonical parameters (3 = 
dS/dU and A = dS/dM z with the microcanonical vari- 
ables (U, M z ) can be obtained by computing the averages 



[3={(3),X = (X),U=(U) and M Z = (M Z ). (107) 



(106) The components of the entropy Hessian 



{kuUi K>um> Kmu-i kmm} can be derived from the gener- 
alized fluctuations relations: 



(Sf3SU) =1 + kuu (SU 2 ) + k um (SM Z SU) , 
(Sf3SM z ) = k uv (SUSM Z ) + n UM (5M 2 ) , 

(6X5U) = kmu (SU 2 ) + kum (6M Z SU) , 
(SX5 M z ) = 1 + kmu (SUSM Z ) + k M m (SM 2 ) . 



(108) 
(109) 
(110) 
(111) 
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by computing the four components of the correlation 
matrix M\ = {(SpSU) , (5/36M z ) , (SXSU) , (SXSM Z )} and 
the four component of the correlation matrix G lJ = 
{(SU 2 ) , (SU5M Z ) , (SMJU) , (SM 2 )}. The microcanon- 
ical entropy S (U, M z ) can be derived from the direct 
numerical integration of the canonical parameters (p, A) 
and the entropy Hessian Kij. The possibility of perform 
such study will be accounted in a forthcoming paper. 

As already shown in the subsection IIV C[ the above 
generalized fluctuation relations can be rewritten by in- 
troducing the generalized response matrix x*- 7 : 



uu 



,MU 



dU 
"dp' 
dM z 
~dp~' 



UM 



MM 



dU 

"ax' 

dM z 



dx ' 



as follows: 



(112) 
(113) 

(114) 



X uu = (SU 2 ) + X UU (S06V) + X UM (SXSU) 
X UM = (5USM) + x uu (6p5M) + x UM (8XSM) , 

(115) 

X MU = (SMSU) + X MU (SPSU) + x MM (SXSU) , (116) 
X MM = (8M 2 ) + X MU (SPSM) + X MM (SXSM) . (117) 



The above expressions exhibit the intrinsic 
reparametrization invariance of the geometrical theory 
developed in the present work. However, conventional 
Thermodynamics usually deals with the total energy 
E = U — BM Z and the external magnetic field B instead 
of the internal energy U and the magnetization M z . 
This is the reason why it is very interesting to rephrased 
the above results in terms of E and B. 

Let us began with the components of the generalized 
response matrix, which can be rephrased as: 



uu 



^ + R^Mi + !/'— + E dMz 
dp dp (3 \dB dB 



M, 



1 fdE dM z , r 



UM 



x MU = 



8M Z 



dp P dB 



1 u dM z 
~xB- 



MM 



1 dM z 



where the symmetric property X — X leads to the 
well-known thermodynamics! identity: 



P dB 



(118) 
(119) 

(120) 



d ^=-p d -Mf-M z =T d -^-M z . (121) 
dB A dp dT K ' 

The fluctuations SU should be rewritten in terms of the 
fluctuations SE, where SU = SE + B5M Z + M Z 5B. The 
ordinary response-fluctuation relations corresponding to 



Sp = SB = are given by: 



<, EE — ^ - /xp2\ -mm _ 1 9M Z /r A/ r2\ 
X ~~dp~ { hX -p~8B- {6Mzh 



-EM 



1 fdE 



X = p{ d -B +M *)=( 5E5M J=X 



-ME 



(122) 
dM z 
~dp~' 

(123) 



where it was introduced what can be considered 
as the components of the response matrix x %1 = 
{x EE ,X EM ,X ME ,X MM } in the representation {E,M Z ). 
The reader can notice that x EE = T 2 C and x MM = 
Txb, where C and xb are the heat capacity and the 
magnetic susceptibility respectively. These classical ex- 
pressions lost their validity when the response functions 
present some anomalous behavior, like C<0orxs<0. 
As already commented, this problem could be solved 
within the fluctuation theory by considering the gener- 
alized response-fluctuation relations when SP ^ and 
SB 0, which can be expressed after some algebra as 
follows: 



x EE 
x EM 



-ME 



-MM 

A 



(SQ 2 ) + X EE (SPSQ) + x EM (SkSQ) , (124) 

(SQSM Z ) + X EE (SPSM Z ) + x EM (SkSM z ) , 

(125) 

(SM Z SQ) + X ME (SPSQ) + X MM (SkSQ) , 

(126) 

(SM 2 ) + X ME (SPSM Z ) + x MM (SkSM z ) , 

(127) 



being SQ = SE + M Z SB and 5k = —PSB, where SQ char- 
acterizes the fluctuations of the transferred heat among 
the interest system and the generalized thermostat. 

The equations (|124j|127j) constitutes the generalized 
response-fluctuation relations in terms of the usual 
control parameters of the conventional Thermodynam- 
ics for a magnetic system, (E, M z \ P, B). The corre- 
sponding generalized expressions for a fluid system is 
very easy to obtain by considering the correspondence 
(E, M Z ;P, B) <-> (E,p; P, V) where p is the pressure and 
V the volume of the container: 



x EE 



x Ep = 



d_E iap 

dp' x pdV 



dE 
dV 



r E = 



dp 
dp' 



(128) 
(129) 



where the symmetric condition x Ep 
well-known thermodynamic identity: 



Y pE leads to the 



X 



dE _ dp 
W ~ ~^dp 



dT 



p. 



(130) 



The generalized response-fluctuation relations are given 
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now by: 

-EE _ I 



X 



)+X EE (6P5Q)+x e p(6kSQ), 
EE (6f35p)+x Ep (SnSp), 
(5pSQ) + x pE (S(3SQ) + X pp {SkSQ) , 
(5 P 2 ) + x pE (5{36p) + x pp (SkS P ) , 

where SQ = SE + p5V and 6k = —j3SV. 



X Ep = (SQSp) 

~pE 

x p 



x pp - 



(131) 
(132) 
(133) 
(134) 



VI. SUMMARY 

It has been considered a generalization of the 
Boltzmann-Gibbs distributions (fTB")) based on the 
reparametrization invariance of the microcanonical en- 
semble. The resulting distribution functions correspond 
to an equilibrium situation where an interest system 
is put in contact with a generalized thermostat in or- 
der to keep fixed the average values of the quantities 

= ( l\ : an external control apparatus becoming 

equivalent for the case of a large enough interest sys- 
tem to the ordinary thermostat with effective fluctuating 
canonical parameters = rj ■ dQ/dl. It is shown that 
the ordinary equilibrium condition between the canoni- 
cal parameters of the thermostat (3 and the correspond- 
ing canonical parameters of interest system (3 — dS/ dl is 
only satisfied now in average (j3^j = //3^ and these ideas 

are used for enhance the possibilities of the well-known 
Metropolis importance sample algorithm in regions char- 
acterized by the ensemble inequivalence. 



The generalized Boltzmann-Gibbs ensemble leads in 
a natural way towards a suitable extension of the clas- 
sical fluctuation theory of the conventional Thermody- 
namics by using a non Riemannian geometric frame- 
work which accounts for the reparametrization changes 
within the microcanonical description. Surprisingly, the 
present approach leads to a novel interpretation of the 
thermodynamic states with negative specific heat and 
other anomalous behavior of the response functions as 
macrostates which can be only controlled by an external 
apparatus with fluctuating control parameters (3. Thus, 
the generalized fluctuation relations act as certain 
kind of thermodynamic uncertainly relations where the 
physical observables of the interest system / and the cor- 
responding effective canonical parameters of the external 
generalized thermostat (3 behave as complementary ther- 
modynamical quantities. The above results constitute the 
basis what could be considered a generalized thermody- 
namic formalism within the microcanonical description. 
A rephrasing the generalized fluctuation relations is car- 
ried out with the introduction of the concept of the con- 
jugate canonical variables (w, I) where uj = (3—0, a repre- 
sentation which possibilities a Riemannian interpretation 
of the fluctuations within the generalized Boltzmann- 
Gibbs description (|12p , constituting in this way a natural 
extension of other geometric interpretations of the clas- 
sical fluctuation theory developed in the past J2^|. As 
example of application, it was obtained the generalized 
response-fluctuation relations for a magnetic system (and 
a fluid system) in terms of the natural control parameters 
of the geometric framework, as well as the usual control 
parameters of the conventional Thermodynamics. 
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